Coherence stabilization of a two-qubit gate by AC fields 
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We consider a CNOT gate operation under the influence of quantum bit-flip noise and demonstrate 
that ac fields can change bit-flip noise into phase noise and thereby improve coherence up to several 
orders of magnitude while the gate operation time remains unchanged. Within a high-frequency 
approximation, both purity and fidelity of the gate operation are studied analytically. The numerical 
treatment with a Bloch-Redfield master equation confirms the analytical results. 
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Despite the remarkable experimental realization of 
qubits Q, Q, El and two-qubit gates 3 in condensed 
matter systems, the construction of a working quantum 
computer remains an elusive goal, not only due to defi- 
ciencies of the control circuitry, but also due to the un- 
avoidable coupling to the environment. Several propos- 
als to overcome the ensuing decoherence have been put 
forward, such as the use of decoherence free subspaces 
[E El Hi IE lp , quantum Zeno subspaces , dynamical 
decouplin g [ill \VA 1 1.1 Il4j , and coherent destruction of 
tunneling 15]. 

A single qubit under the influence of bit-flip noise 
can be modeled by the spin-boson Hamiltonian Hq = 
— ye 2 + cr x £, where a x,z denotes Pauli matrices and £ 
is a shorthand notation for the quantum noise specified 
below. The influence of the bath is governed by the 
spectral density of the noise at the tunneling frequency 
A/h. A possible driving field may couple to any pro- 
jection n of the (pseudo) spin operator a, i.e., be pro- 
portional to n ■ a. In Ref. [H], two particular choices 
have been studied and compared against each other: A 
driving of the form H (t) = Aa z cos(r2i) commutes with 
the static qubit Hamiltonian while it modifies the bath 
coupling a x t; in such a way that the spectral density of 
the bath at multiples of the driving frequency becomes 
relevant. For a proper driving amplitude, this eliminates 
noise with frequencies below the driving frequency which, 
thus, should lie above the cutoff frequency of the bath. 
This scheme represents a continuous- wave version of dy- 
namical decoupling [lflj . By contrast, a driving of the 
type H(t) — Ao x cos(Oi) renders the qubit-bath cou- 
pling unchanged but renormalizes the tunnel splitting A 
towards smaller values and thereby causes the so-called 
coherent destruction of tunneling (CDT) 0, . Then, 
decoherence is determined by the spectral density of the 
bath at a lower effective tunnel frequency. For an ohmic 
bath being linear in the frequency, the consequence is 
that both decoherence and the coherent oscillations in 
the rotating frame are slowed down by the same factor 
[lfij . Therefore, the number of coherent oscillations is 
not enlarged and, thus, for single-qubit operations, CDT 
might be of limited use. 

In this work, we propose a coherence stabilization 



scheme for a CNOT gate based on isotropic Heisenberg 
interaction 0, llfi| . Our scheme does not suffer from 
the drawbacks mentioned above because (i) it involves 
only intermediately large driving frequencies that can lie 
well below the bath cutoff and (ii) it does not increase 
the operation time of the gate. Since we shall employ a 
field that couples to the same coordinate as the quantum 
noise, the present coherence stabilization is different from 
dynamical decoupling. 

CNOT gate with bit- flip noise. — We consider a pair of 
qubits described by the Hamiltonian |l£|,L!j4|20, 21] 
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with a spin-spin coupling of the Heisenberg type, where 
j labels the qubits. In order to construct a quantum 
gate, the tunnel splittings Aj, the biases 6j, and the 
spin-spin coupling J have to be controllable in the sense 
that they can be turned off and that their signs can be 
changed. Then, a suitable se quen ce of interactions yields 
the CNOT operation [3 H Hi El 



U, 



H 

CNOT 



Uh(tt/8)U 1z (tt)Uh{tt/8). 



(2) 



where Ui z (tp) = exp(—i(po~ z /2) represents a rotation of 
qubit 1 around the z-axis and Uh{^) = exp(— itpai ■ a%) 
describes the time evolution due to the qubit-qubit in- 
teraction. The symbol ~ denotes equality up to local 
unitary transformations, i.e., transformations which act 
on only one qubit. Since we focus on decoherence dur- 
ing the stage of the qubit-qubit interaction, we take as 
a working hypothesis that one-qubit operations perform 
ideally Thus, decoherence takes place while Aj = 6j = 
and J > during the required total qubit-qubit interac- 
tion time tj — itK/AJ. 

The bit-flip noise is specified by the system-bath 
Hamiltonian H = H gate + H coup i + iibath where H coup i = 

\ Si=i,2 a j ft c "( a ]u + a ]v) describes the coupling of 
qubit j to a bath of harmonic oscillators with frequen- 
cies w„ and /?bath = Y^rjv h^v&\ v a jv The bath is fully 
specified by its spectral density I(u>) — c^5{uj — uj v ). 
Within the present work, we consider the so-called ohmic 
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spectral density I(oj) — 2-7rau;exp(— lo/uj c ) with the di- 
mensionless coupling strength a and the cutoff frequency 
lu c . In order to complete the model, we specify the 
initial condition for the total density matrix to be of 
the Feynman- Vernon type, i.e., initially, the bath is 
in thermal equilibrium and uncorrelated with the sys- 
tem, Ptot(*o) = p(to) ® -Rbath,oq, where p is the re- 
duced density operator of the two qubits and -Rbath.cq oc 
cxp(— iJbath/fcsT) is the canonical ensemble of the bath. 

If aksT is larger than the typical system energy J 
and the dissipation strength is sufficiently small, a -C 1, 
the dissipative system dynamics is well described within 
a Born-Markov approach. There, one starts from the 
Liouville-von Neumann equation ihp tot = [H, p tot ] for 
the total density operator and obtains by standard tech- 
niques the master equation 

P = |[^ate, p] - YMl Wi(*)> PW - EK' i P M P}\ 

3 3 

= ^[H gate ,p]-A(t)p (3) 
with the anti-commutator {A, B} = AB + BA and 

Qj(t) = T- / dT / duS(oj)cos(uJT)a?(t-T,t). (4) 
4tt J q J 

Here, <S(w) = I{uj) coth(fiw/2fcsT) is the Fourier trans- 
formed of the symmetrically-ordered equilibrium bath 
correlation function §({£(t), £(0)}} cq of the collective co- 
ordinate £ = c„(aj, + a„). The notation X(t, t') is a 
shorthand for the Heisenberg operator W (t,t')XU(t,t') 
with U being the propagator of the coherent system dy- 
namics. Note that S(u) and I(lu) are the same for both 
qubits due to the assumption of two identical environ- 
ments. Replacing in Eq. (@J the term S(u) cos(wt) by 
I (to) sin(wT), yields the operator Pj(t). We emphasize 
that the particular form J21 °f the master equation is 
valid also for an explicitly time-dependent Hamiltonian. 

The heat baths, whose influence is described by the 
second and third term of the master equation J3J), lead 
to decoherence, i.e., the evolution from a pure state to 
an incoherent mixture. Decoherence can be measured by 
the decay of the purity tr(p 2 ) from the ideal value 1. The 
gate purity (later refered to as "purity" )V(t) = tr(p 2 (i)), 
which characterizes the gate independently of the specific 
input, results from the ensemble average over all pure ini- 
tial states 23] . For weak dissipation, the purity is deter- 
mined by its time derivative at initial time. Performing 
the ensemble average, one finds the gate purity decay 

V(t)\ t=Q = -2trW= --^JL^TrA, (5) 

where d — 4 is the dimension of the system Hilbert space 
and Tr denotes the trace in superoperator space defined 
as Tr A = E a ,a' tr (l a ')( a l A l a )( a 'l) with {»} bein S an 



arbitrary orthonormal basis of the system Hilbert space. 
It can be shown that the term containing Pj does not 
contribute to Tr A. 

In order to evaluate for our setup the purity decay, we 
need explicit expressions for the operators Qj and, thus, 
have to compute the Heisenberg operators ait — r, t) for 
the Hamiltonian Hq = Ja\ ■ a-i . This calculation is most 
conveniently done in the eigenbasis of Hq, i.e. in the basis 
of the total (pseudo) spin L = \{(Ti +^2). After perform- 
ing the time and the frequency integration in Eq. we 
finally arrive at V = -|{5(0)+5(4J/fi)} 5 where we have 
ignored Lamb shifts and defined 5(0) = lim^^o S(lu) = 
inakeT '/h. In particular, we find that for low temper- 
atures, fcflT < J, decoherence is dominated by S{AJ/K) 
such that V w —W%aJ /5h. This value reflects the in- 
fluence of the so-called quantum noise which cannot be 
reduced by further cooling. 

AC driving field. — In order to manipulate the coher- 
ence properties, we act upon qubit I by an ac field which 
causes a time-dependent level splitting according to 

H ac (t) = f(t)af, (6) 

where f(t) is a T-periodic function with zero time- 
average. Due to the finite driving time, which acts only 
while the Heisenberg coupling is switched on, the spec- 
trum of the driving field acquires a dispersion Aw w Q/k, 
where k is the number of driving periods. To keep the in- 
fluence of the dispersion small, we have to choose k 3> 1. 

Next, we derive within a high-frequency approximation 
analytical expressions for both the coherent propagator 
U(t,t') and the purity decay 10. We start out by trans- 
forming the total Hamiltonian into a rotating frame with 
respect to the driving by the unitary transformation 

U Be (t) = e- i *W<*, m = \f dt' f{t'), (7) 

which yields the likewise T-periodic Hamiltonian H (t) = 
Ul c (t)H gatc U ac (t) . For large driving frequencies f2 3> J, 
it is possible to separate time scales and thereby replace 
H(t) by its time-average 

H = (J - JxKof + Jx ff x ■ (8) 

where the constant J± — J(cos[2<f>(t)])T denotes an ef- 
fective interaction "transverse" to the driving and (■ ■ -)t 
the time average over the driving period. Within this ap- 
proximation, the propagator of the driven system reads 

U eS (t, t') = e -*C*K e -iH(t-t')/h e #(t'K . (9) 

Having this propagator at hand, we are in the position to 
derive explicit expressions for the operators crj (t — t, t) 
and Qj and, consequently, also for the generator of the 
dissipative dynamics A. The calculation is conveniently 
done in the basis of the total spin L and L x which, owing 
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to the relation erf erf = \ (erf +erf ) 2 — 1, is an eigenbasis of 
the Hamiltonian (JSJ. We insert the resulting expression 
for A into Eq. © and finally obtain the purity decay 

V = -\{S{Q)+S{Ujh)}. (10) 

For /(f) = 0, wc find Ji = J such that Eq. (fTUf) agrees 
with what we found in the static case; otherwise, the 
inequality | Jj_| < J holds and, thus, the bath correlation 
function S in Eq. IjlOjl has to be evaluated at a lower 
frequency. For an ohmic or a super-ohmic bath, S{ui) 
is a monotonously increasing function and, consequently, 
the ac field reduces V (unless J > w cuto ff)- 

The purity decay assumes its minimum for J± = 0. 
This condition marks the working points on which we 
shall focus in the following. For an ohmic spectral density 
I(lo) = 2nauj, the purity decay at the working points be- 
comes V — — !<S(0) = —8irak B T/5 for all temperatures. 
This particular value has to be compared to the purity 
decay in the absence of driving: An analysis reveals that 
for k B T > J, decoherence is essentially driving indepen- 
dent. By contrast for low temperatures, fcgT < J, the 
driving reduces the decoherence rate by a factor k B T/2J. 
The reason for this low-temperature behavior is that for 
J± = 0, the coupling operators erj commute with the ef- 
fective Hamiltonian (JBJ) such that noise acts as pure phase 
noise with a strength proportional to the temperature. 

For a rectangular driving for which /(f) switches be- 
tween the values ±A/2, the condition J± = yields 
A = Ml which corresponds to the application of two 
7r-pulses per period. For a harmonic driving, /(f) = 
Acos(Of)/2, one obtains J± = JJo(A/t0), where Jo de- 
notes the zeroth-order Bessel function of the first kind. 
Then, at the working points J± = 0, the ratio A/ Ml as- 
sumes a zero of Jo, i.e., one of the values 2.405.., 5.520.., 
8.654.., .... 

So far, we focussed on decoherence and ignored the in- 
fluence of the driving on the coherent dynamics. Since 
the driving changes also the coherent dynamics, the pulse 
sequence of the CNOT operation needs a modification: 
At the working points of the driven system, the prop- 
agator © becomes U e g(t,t') = exp[— iJof erf (t — t')/H] 
and, thus, represents the time evolution caused by a so- 
called Ising interaction Jafa%. This time evolution al- 
lows one to implement the alternative CNOT operation 
^cnot ~ exp(-i7rcrfcrf/4) = U cS (t + tj,t) HJUj. Note 
that the interaction time tj = irh/AJ is the same as for 
the original gate operation Uq NOT . Since U ac (2ir/fl) is 
the identity [cf. Eq. 10] , we assume for convenience that 
the driving period 27r/f2 is an integer multiple of f j, i.e. 
fl = 8kJ/H with integer k . 

Numerical solution. — The calculation of the purity 
decay presented above relies on a high-frequency ap- 
proximation which is correct only to lowest order in 
J /Ml. Thus, these should be compared to the exact 
numerical solution of the master equation @. An ef- 
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FIG. 1: Purity loss 1 — V during the Heisenberg interaction 
time tj as a function of the driving amplitude for various 
temperatures. The driving frequency is f2 — 32J/H and the 
dissipation strength 2-na — 0.01. For A = 0, the undriven 
situation is reproduced. 



ficient scheme for that purpose is a modified Bloch- 
Redfield formalism whose cornerstone is a decomposi- 
tion into the Floquet basis of the driven system [22| : 
According to the Floquet theorem, the Schrodinger 
equation of a periodically driven quantum system pos- 
sesses a complete set of solutions of the form [^(f)) = 
exp(— ie a t/h)\<p a (t)}, where the so-called Floquet states 
\4> a (t)) obey the time-periodicity of the Hamiltonian and 
e Q denotes the so-called quasienergy. They are elements 
of an Hilbert space extended by a T-periodic time coor- 
dinate and are computed from the eigenvalue equation 
[H(t) - id/dt]\(j) a (t)) = e a \(f> a (t)}. In the Floquet ba- 
sis {\4>a(t))}, the master equation © obtains the form 
Pap = — - ep)Pap ~ E a '/3' ^-ap.a'P' Pa' p> ■ The ben- 
efit of this representation is that the time dependence of 
the original Hamiltonian has been eliminated by choos- 
ing a suitable basis. Moreover, for weak dissipation, we 
can replace within a rotating-wave approximation A(f) 
by its time average j2^ ■ Finally, we integrate the master 
equation to obtain the dissipative propagator which in 
turn allows to evaluate all quantities of interest. 

In our numerical studies, we restrict ourselves to purely 
harmonic driving /(f) = Acos(57f)/2. The resulting pu- 
rity loss during the interaction time tj is depicted in 
Fig. ^ We find that for k B T > J, decoherence is fairly 
independent of the driving. This behavior changes as 
the temperature is lowered: Once k B T < J, the purity 
loss is significantly reduced whenever the ratio A/Ml is 
close to a zero of the Bessel function Jo. Both observa- 
tions confirm the preceding analytical estimates. The be- 
havior at the first working point A w 2 A Ml is depicted 
in Fig. Ufa). For relatively low driving frequencies, we 
find the purity loss being proportional to J/Cl. There 
are significant deviations from the analytical result for 
small because the low-frequency regime is not within 
the scope of our analytical treatment. With increasing 



4 



p. 

i 



i 



10" 



10" 



10" 



10" 



10" 



10" 



I ' — I ' — I ' — I 



oc j/n 



— k B T = 0.1J 
---- k B T = 0.01J 

. /••;;/ = 0.001J ...~.'.~'.--.r 

I I 1 1 ll| 1 — I I I 1 1 ll| 1 — I I I I ll| 1 — I I I 1 1 ll| 1 — ! 



(b) 



1_ 



_l I ' ■ ■ ' 



10 1 10 2 



10 d 



10 4 



n [j/h] 



FIG. 2: (a) Purity loss during the Heisenberg interaction as 
a function of the driving frequency. The driving amplitude 
is adjusted such that 1 — V assumes its first minimum; i.e. 
A ~ 2.4 HQ,, cf. Fig. The dotted lines mark the analytical 
estimate 1 — V(tj) w ~V(0)tj and the proportionality to 
J/Q. (b) Corresponding fidelity defect 1 — T. 



driving frequency, the discrepancy decreases until finally 
thermal noise dominates decohcrence. The numerical so- 
lution confirms almost perfectly the results from the high- 
frequency approximation. 

Still, there remains one caveat: The gate operation 
^CNOT relies on the fact that U e g is a good approxima- 
tion for the dynamics of the driven gate because any devi- 
ation represents a coherent error. Therefore, we still have 
to demonstrate that such coherent errors are sufficiently 
small. As a measure for t his, we employ the so-called fi- 



delity T = tr[pideai p(tj)] [23| which constitutes the over- 
lap between the real outcome of the operation, p(tj), and 
the desired final state pideai = U I (Tr/A)p in uj(TT/A) in the 
average over all pure initial states. An ideal operation is 
characterized by T = 1. Here, Ui(ip) = exp(— iipaf ) 
denotes the propagator due to an ideal Ising spin-spin 
interaction. Figure H£b) shows the fidelity defect 1 — T 
of the real gate operation at the first working point. In 
the relevant temperature regime k B T < J, we find the 
crucial result that the fidelity defect is of the same or- 
der or even smaller than the purity loss. Thus, we can 
conclude that coherent errors are not of a hindrance. 

In summary, we have put forward a strategy for re- 
ducing decoherence in a two-qubit quantum gate during 
the stage of qubit-qubit interaction. This coherence sta- 
bilization scheme is different from previous proposals in 
two respects: First, it does not rely on the manipulation 



of the coordinate that couples the qubits to the envi- 
ronment. By contrast, the central idea of our scheme is 
rather to suppress the coherent system dynamics "trans- 
verse" to this sensitive system coordinate. In that way, 
the bit-flip noise, which produces unwanted transitions, is 
turned into pure phase noise which is proportional to the 
temperature. This enables a coherence gain by cooling. 
The second difference is that the proposed scheme elimi- 
nates also the noise coming from the spectral range above 
the driving frequency and, thus, is particularly suited for 
ohmic noise spectra with large cutoff frequencies. More- 
over, the driven system still allows one to perform the 
desired CNOT operation with high fidelity and within 
the same operation time as in the absence of the control 
field. Hence, the gained coherence time fully contributes 
to the number of feasible gate operations. 

This work has been supported by the Freistaat Bayern 
via the quantum information initiative "Quanteninfor- 
mation langs der A8" and by the Deutsche Forschungs- 
gemeinschaft through SFB 631. 
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